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Abstract 

The theme of this paper is the applications of generalized function (GF) theory to the wave equation in 
aeroacoustics. We start with a tutorial on GFs with particular emphasis on viewing functions as continu- 
ous linear functionals. We next define operations on GFs. The operation of interest to us in this paper is 
generalized differentiation. We give many applications of generalized differentiation, particularly for the 
wave equation. We discuss the use of GFs in finding Green’s function and some subtleties that only GF 
theory can clarify without ambiguities. We show how the knowledge of the Green’s function of an 
operator L in a given domain I) can allow us to solve a whole range of problems with operator L for 
domains situated within D by the imbedding method. We will show how we can use the imbedding 
method to find the Kirchhoff formulas for stationary and moving surfaces with ease and elegance with- 
out the use of the four-dimensional Green’s theorem, which is commonly done. Other subjects covered 
are why the derivatives in conservation laws should be viewed as generalized derivatives and what are 
the consequences of doing this. In particular we show how we can imbed a problem in a larger domain 
for the identical differential equation for which the Green’s function is known. The primary purpose of 
this paper is to convince the readers that GF theory is absolutely essential in aeroacoustics because of its 
powerful operational properties. Furthermore, learning the subject and using it can be fun. 

Nomenclature 


C°° : The space of infinitely differentiable functions 

C n : The space of continuously differentiable functions up to the order n 
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D : The space of infinitely differentiable functions with compact (closed and bounded) support 
D' : The space of continuous linear functionals on the space D, i.e., the space of Schwartz GFs 
GD : Generalized derivative 
GF : Generalized function 

OF : Ordinary function, any function that is locally Lebesgue integrable 
e : Belongs to 

[0, 1] \ [a, b] : Set subtraction, [0, a] U ( b , 1] 


1- Introduction 


The linear wave equation and its counterpart in frequency domain, the Helmholtz equation, are the most 
important partial differential equations of aeroacoustics. For example, both the Lighthill 1 jet noise 
equation and the Ffowcs Williams-Hawkings (FW-H) equation 2 are linear wave equations. We will 
primarily be concerned with the wave equation here. Although this equation in one, two, three and 
higher space dimensions has been analyzed since the eighteenth century 3 , there are still many issues and 
subtleties associated with finding the solution, particularly in aeroacoustics, that are in need of further 
attention. It is true to say, and this fact is well-known to mathematicians, that the appropriate mathemati- 
cal tool for solving PDEs is the theory of distributions or generalized functions (GFs) 3-12 . The deriva- 
tion and finding the solution of the FW-H equation depend heavily on the concept of GFs which 
requires much mathematical maturity from the users. Furthermore, even when we use the Lighthill 
equation, we may be dealing with flowfields that have discontinuities within them, such as shocks, 
which can be sources of sound. The natural tool to work with these discontinuities is GF theory. The 
principal aim of this paper is to give a working knowledge of the multidimensional GF theory to acousti- 
cians, particularly when using the linear wave equation. 

There are many books and papers on the mathematics of the wave equation and the related problems 
that are used by acousticians. Unfortunately, almost all of the commonly referenced books and papers 
were written before, or without the use of, the theory of distributions (generalized functions) of Laurent 
Schwartz 4 . These books and papers use an intuitive approach to Dirac delta function and its algebra that 
is not rigorous and is prone to giving erroneous results. Furthermore, much of the true power of GF 
theory is not apparent from the intuitive approach. The available books on mathematics that discuss 
generalized functions (GFs) rigorously 4- 6 are also beyond the reach of most acousticians because of 
the level of abstraction involved to comprehend the theory. We believe that most of the abstractions are 
not necessary for application of the GFs. The situation is much like the use of the real numbers by the 
man in the street in everyday business without ever knowing about Cauchy convergent sequences 
because he only needs to know the useful properties of the real numbers. For mathematically minded 
persons the abstractions in the development of the theory can be learned most easily after a working 
knowledge of GFs is attained. 

We will start by giving the definition of GFs as continuous linear functionals on a given space of func- 
tions. We will do this in the simplest possible way. We will then discuss how the definition of many 
useful operations on ordinary functions can be extended to GFs. Of particular interest to us here is the 
definition of generalized derivative and its many useful applications. These include the Green’s func- 
tions of ODEs and PDEs, finite part of divergent integrals, the Kirchhoff formula for moving surfaces 
and getting the jump condition from conservation laws. We will discuss some questions associated with 
the derivation of the FW-H equation. In particular, we answer the question of whether any function 
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other than the null function can be assumed inside the moving surface in the derivation of the FW-H 
equation. We will see that in every aspect of working with the linear wave equation, the natural mathe- 
matical tool is the GF theory. 

We wish to convey to the readers our enthusiasm of working with GFs as the most essential, useful and 
powerful tool for application in acoustics. We have, therefore, chosen a subtitle for this paper with an 
approving nod to the title of the delightful and scholarly book of Thomas Korner 17 . 

2- Generalized Functions (GFs) 

In this section we will first start by giving the definition of GFs. We then present the most common 
operations on GFs. We will pay a special attention to generalized differentiation which is the fundamen- 
tal concept for obtaining the Green’s function, finding the finite part of divergent integrals and working 
with jump conditions from conservation laws. We will then give many applications of GFs relevant to 
acoustics. We will say little about generalized Fourier transformation because that subject by itself 
requires a full paper. 

In learning GF theory, you should know that there are several different approaches to the definition of 
GFs. The method we present here is the functional approach of Schwartz 4-11 . The sequential 
approach of Temple and Lighthill 14 , Jones 1 2 and Mikusinski 15 . The completeness theorem of GFs 
relates the two approaches 7 . The approach by Bremermann defines GFs as the boundary value of ana- 
lytic functions 16 . This subject is valuable in studying generalized Fourier analysis and, in its full general- 
ity, requires the knowledge of the highly specialized topic of functions of several complex variables. 
Most of Bremermann’ s book 16 does not require this topic and for learning about one dimensional GFs, 
it is easily accessible to readers with a good knowledge of analytic functions of a single complex vari- 
able. Eventually, when you work on advanced Fourier transform theory, you will have to learn 
Bremermann’s approach because it sheds light on the derivation of many of the results obtained by 
physicists and applied mathematicians using ad hoc and mostly intuitive reasoning. The subject can be 
shown to be also closely related to the functional approach we use here but we will not say any more 
about this connection 7 . 

The Polish mathematician Jan Mikusinski introduced one more method of defining GFs known as the 
operational approach which is based on the concepts from abstract algebra 17,18 . Mikusinski’s publica- 
tions are beautifully written and are easy to read and comprehend. Apart from explaining rigorously the 
Heaviside operational calculus, this approach can also be used to solve an entirely different set of prob- 
lems than the first three approaches, such as finding the exact solutions to finite difference equations and 
recursion relations. The relationship between the Mikusinski’s operational approach and the functional 
approach cannot be established easily and we are not aware of any publications on the subject. We 
encourage the readers to study this approach because of its usefulness and power in solving some impor- 
tant problems of physics and engineering 17- 20 . We will say no more on this approach here. 

And finally, the readers should be aware of the fact that in recent years GF theory has been extended to 
the so-called nonlinear GF (NLGF) theory where some of the limitations of Schwartz distribution (GF) 
theory are removed 21- 22 . For example, it can be shown that one cannot extend the operation of multipli- 
cation to all of Schwartz GFs. This shortcoming of the Schwartz theory does cause some difficulties in 
applications that Colombeau 21 , Rosinger 22 and Oberguggenberger 27 as well as other mathematicians 
have tried to resolve. Currently, the NLGF theory is very abstract and difficult to master. The subject is 
under intense development but is badly in need of simplification. However, it appears that the subject of 
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nonstandard analysis (NSA) 24 will eventually be able to clarify and simplify NLGF theory considerably 
and bring the theory within the reach of physicists and engineers 25 . We will not cover this subject here. 
The readers should know that to understand NLGF theory and NSA, requiring one to have a knowledge 
of several esoteric areas of mathematics including topology, abstract algebra and mathematical logic. 
Fortunately, there are many expository books and articles on these subjects on the internet, particularly 
on NSA. We will not cover this theory here. 

For those interested in the history of GF theory, we recommend the books by Liitzen 26 and Dieudon- 
ne 27 as well as the article by Synowiec 28 . In the opinion of Liitzen, the Russian mathematician Sobolev 
invented the concept of GFs (more exactly, generalized derivatives) and Schwartz developed the theory 
of GFs without being aware of Sobolev’s earlier contributions 4 . As expected, many other mathemati- 
cians contributed to the field. Without a doubt, the members of the Moscow School of Mathematics 
under I. M. Gelfand have contributed substantially to the development and applications of GF theory 
making the subject accessible to physicists and engineers. The well-known book on GFs by Gelfand and 
Shilov 7 is the first of a six volume opus by this School published in the West in the sixties of the last 
century. This monumental work, which is still read today, has popularized the GF theory enormously 
among physicists and mathematicians by its breadth and depth as well as by the clarity of the exposition. 
For learning about more recent advances we suggest the books by Hormander 5 and Taylor 6 which are 
both written for mathematicians. 

Most of this section is based on the first author’s expository NASA publication 29 . In addition to this 
publication, we recommend that you consult the books by Strichartz 10 , Stakgold 8 , Gelfand and Shilov 7 , 
and Kanwal 11 . For beginners, we particularly recommend the excellent book by Strichartz for deeper 
understanding of many of the concepts discussed here. This book emphasizes GFs in one dimension 
which is essential in some applications, e.g., time series analysis. However, there are many problems of 
physics and engineering for which multidimensional GF theory is very important 7-9,11,12,29 . This is 
the subject that we concentrate in the present paper. The emphasis will be on using the theory to solve 
problems rather than stating theorems and giving proofs. 

2.1- Definition of Generalized Functions (GFs) 

What we call a function, say fix) , in calculus can be thought of as a table of ordered pairs (x, fix)). We 
can visualize this table as a graph of x versus y = fix) . If we want to generalize the concept of a func- 
tion, we must change the definition of a function in such a way that our old functions are still considered 
functions and we have many new useful mathematical objects that are functions by the new definition. 
Once we introduce a new definition of a function, we must make sure that we can extend all the useful 
operational properties of the old functions, e.g., addition of two functions, to all the new mathematical 
objects which we now collectively (the old functions and the new objects) call generalized functions 
(GFs). For now, we refer to all the functions that we learned about in calculus as ordinary functions 
(OFs). We will elaborate further on this below. 

Let us start by considering an ordinary function g(x) which is periodic over the interval [0, 2n\. We 
know that the complex Fourier components of this function are given by the relation 

l r 2n 

G n = I g (x)e inx dx (w = — oo to oo) 

Jo 


( 1 ) 



FF, KSB, FF-AIAA-CEAS 2007F.nb 


5 


We know that g(x) = E”=-oo G n e inx , that is, the periodic function g(x) can be constructed from the 
knowledge of the table with infinite entry {G„, n - -ootooo}. From another point of view, the set of 
complex numbers {G„} can be considered as the mapping of the space of functions [ e' n x , n = -oo to 00} 
into complex numbers by the rule given in eq. (1). Such a mapping from a given space of functions into 
scalars (real or complex numbers) is called a functional. Remember that a functional is always a rule 
that describes the mapping from a given space of functions into scalars. The numbers { G n } are the 
functional values of this functional. Note that the function g(x) enters the definition of the functional. As 
a matter of fact, we should use the symbol G n [g\ to remind us that the number G n depends on the func- 
tion g{x) . So far we have shown that the periodic function g(x) can be described (identified, defined) by 
the table of its functional values: 



( 2 ) 


We have assumed here that the function g(x) has a Fourier series (e.g., satisfies the Dirichlet conditions, 
etc.). Now we have the germ of an idea to extend the definition of a function. 

Let us take a space of functions of one variable D with a large number of well-defined functions in the 
space. We will say more about the properties of the functions in this space soon. Let us select an ordi- 
nary function f(x) and for any function <p (x) in the space D, define the functional (i.e., the rule): 

/>oo 

F[(f>] = I f{x)(p{x)dx (/ is fixed in this rule and f varies all over D) (3) 

— OO 

Note carefully how the function fix) enters the definition of the functional. To describe /, this function 
is kept fixed in the rule given by eq. (3) and different functions (f> are taken ( conceptually , of course, not 
literally ) from the space D to form the table \F\(/)\, (f> e D\ . Therefore, we must have so many functions 
in space D that any two different ordinary functions generate different tables of functional values at 
least for some entries. The space D is known as the test function space. The question is how to select 
the space D in such a way that the set of the functional values {F[<p\, f sD} represents (i.e., describes or 
has as much information content as) any ordinary function f(x) . Before we do this, we must now elabo- 
rate on what we mean by ordinary functions. 

By ordinary functions (OFs), we mean the locally Lebesgue integr able functions. This means that they 
are functions that are Lebesgue integrable 7- 11,19 over any finite interval. These functions include any 
conceivable functions, no matter how wild, that are needed in applications. If you do not know Leb- 
esgue integration, you can think of Riemann integration that you learned in calculus courses. Lebesgue 
integration is a generalization of the Riemann integration that can be used for the integration of most 
functions that do not have Riemann integral. 

One can show that there is an uncountable number of OFs. Therefore, the space of test functions must 
contain an uncountable number of functions so that any two different OFs can be distinguished from 
each other when viewed as tables of functional values. This means that our tables of functional values 
used to represent OFs are also uncountable. Furthermore, since OFs are as wild as they can be, we 
require that the functions in space D to be as well-behaved as they can be so that the integral in eq. (3) is 
always defined. This leads us to require that the functions in D to be infinitely differentiable (denoted 
C°° functions). Also to ensure the integral in eq. (3) does not become infinite, we assume that all func- 
tions in space D are identically zero outside of a finite interval. This interval can be different for differ- 
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ent test functions. A function that is identically zero outside some finite interval is said to have bounded 
support. Thus, we have come to the following conclusion regarding the space of test functions D that we 
need to use to identify OFs by their tables of functional values: 

The test function space D that one needs to use in defining OFs by their tables of functional 
values is the space of all C°° functions with compact ( closed and bounded) support. 

Let us say something about this space next. The best known, and in many ways the fundamental, exam- 
ple of a function in D is the following bell-shaped function with the peak value at the origin. It is 
defined for any real number a > 0 : 


0 (a; a) 


( g x 2 /(x 2 -a 2 ) 

\x\ < a 

1° 

| x | 5= a 


( 4 ) 



“t 

-a 


Figure 1- The Graph of the function 4> (jc; a) 

It is easy to see that the support of this function is [- a, a] and it is infinitely differentiable on the open 
interval (- a, a ) . Figure 1 shows the graph of this function. The proof of infinite differentiability at 
points x = - a and x = a is much harder and is an exercise in calculus albeit requiring much technical 
skill. You should know that at these two points all the derivatives of this function are zero. Are there any 
other types of functions in space D besides d> (x; a) ? The answer is affirmative and we will turn to this 
point next. 

One proves rigorously in real analysis that there is an uncountably infinite number of continuous func- 
tions. For any continuous function g(x) defined on the bounded interval I = [b, c] , the convolution of 
g(x) and ([> (x; a) on I gives a C°° function with the bounded support [b - a, c + a) 19 . Therefore, there 
is an uncountable number of functions in the test function space D as we are required to have. Another 
fact to remember about this space is that the functions in space D are not analytic even though they are 
C°° functions! 
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At this stage, we have a new way of looking at OFs but our original intention was to extend the defini- 
tion of a function so that we have new objects that are not OFs but are considered functions. To do this 
we study two important properties of the functionals described by eq. (3). The first one is very easy to 
see: 

1 - The functional F\d> \ is linear. 

This simply means that for a and / 3 scalars and 0 and ip functions in D , we have 

F[a f + j3 if] = a F[f\ + /3 F[if] (5) 


To show this, replace d) (x) on the right of eq. (3) with ot(p + fiif and use the linearity property of inte- 
grals followed by functional interpretation of the result to get the right side of eq. (4). We remind the 
readers that the OF f(x) used in the definition of the functional of eq. (3) is fixed here and the functional 
F\(/>\ is, in point of fact, representing f(x) . 

The second property is much harder to explain and comprehend. It is: 

2 - The functional F[0] is continuous. 

The first thing you must know is that this property, whatever it means, gives GFs a lot of useful proper- 
ties that we need in applications. Second, to the beginners the definition of continuity of a functional 
appears too complicated and perhaps a futile attempt in abstraction. We admit that the definition of the 
continuity of a functional is complicated but its definition cannot be stated in simpler form than we give 
here. The original definition by Schwartz was much more complicated and based on locally convex 
topological spaces 4 . Just remember that you will appreciate the full usefulness of this property when 
you learn some more about GFs and their applications. See, for instance, Example 3 below. 

We first distinguish some sequences of functions {0„}in D, which we say they go to zero in D, and 

D 

write it as (p n — > 0. For any such sequence, we require that it satisfies two conditions: 

i- There is a bounded interval J such that the support of (p n c J for all n, and 

ii- (p n and all its derivatives go to zero uniformly in J . 


We make two useful comments here concerning these conditions. First, in general, different sequences 
going to zero in D have different bounded interval J associated with them. In other words, we require 
only that for a given sequence going to zero in D such a bounded interval to exist. Second, condition ii 
simply means that if m is the order of the derivative, m = 0, 1 , 2, ...., and given 
s > 0, we can find N(m) such that for any x in the interval J, we have \(p^ (x)| < e for all n > N{m). 
Note that the word “uniformly” means that the number N(m) is independent of x but it can depend on 
the order of differentiation m . Here we have used the notation (pff (x) for the mth derivative of (p n (x) . 

Example 1- We can show easily that the following sequence goes to zero in D : 

(p n (x) = — 0 (a; a) (See eq. (4) for the definition of 0 (x; a)) (6) 

n 


Example 2- The following sequence of functions does not go to zero in D : 

1 ( x \ 

c p n (x ) = — 0 — ; a) (See eq. (4) for the definition of 0 (x; a)) 
n V n ' 


(7) 
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The reason is that for a given function (f> n , the support is the interval [-na, na] and, therefore, the 
support cannot remain bounded for all n. Condition ii is, thus, violated. 

Now we can give the definition of continuity of a functional: 

D 

A functional F[cf)] is continuous if for any d> n — > 0, we have F[d> n ] — * 0 . 

Note that { F\(f> n | } is simply a sequence of real or complex numbers. By using the properties of (Leb- 
esgue) integrals, we can show the following very important result: 

Any OF fix) generates a continuous linear functional through 

X oo 

f(x) 4>{x)dx (d> is in the test function space D) 

“OO 




F[<W=/f(x) 4>(x) dx 

Figure 2- The functional rule used to define an OF f{x) by a table of its functional values over the 
space D 


Figure 2 shows schematically how the table of functional values is constructed for an OF. So we can 
use the above functional (or rule) to construct the sequence or the table of numbers {F[<A], d> e D\ to 
identify (or describe , distinguish) the OF fix) . Are there any more continuous linear functionals 
described by other rules? The answer is affirmative and, in fact, there are many of them. Let us consider 
the following example: 

Example 3- Define the functional 5 \d)\ = <l> (0), feD. We can easily show that this functional is linear. 
We now show that it is also continuous. Take any sequence that goes to zero in D, say {f n }. Since 

D 

f n — > 0, then in particular ([>,, — >0 uniformly. The uniformity implies that it does not matter where x is. 
Therefore, (0) = 6 \ d>„ | — >0. Can you see how we use condition ii of the definition of continuity of a 
functional here? Because of the importance of this functional in GF theory, we will give it the name of 
the Dirac functional. This is obviously related to the Dirac delta function. We will establish the connec- 
tion below (see eq. (9)). 

Using elementary Lebesgue integration theory, we can show that there is no OF 6(x) that can generate 
the continuous linear functional 6 \d)\ by using the rule given in eq. (3). Therefore not all the continuous 
linear functionals are generated by OFs based on eq. (3)! 

Now if we decide to define (or describe, identify) any OF conceptually by the table of its functional 
values using the rule of eq. (3), we may just as well say that the table generated by any continuous 
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linear functional also defines a function. We now have new objects (tables) that cannot be identified 
with OFs, e.g., the functional in Example 3. We call the entire collection of objects (tables of functional 
values) generalized functions (GFs). We have achieved our goal of generalizing the definition of OFs. 
We, thus, state the following: 

Generalized functions are defined by the tables of functional values of continuous linear functionals. 

Since we are now going to think of functions in the new way by their functional values, we will use a 
slight abuse of language and say: 

Generalized functions are continuous linear functionals. 

The particular space of GFs that we have constructed using the test function space D is known as the 
(Schwartz) GF space D ' . We know that by construction, all OFs are also GFs. We refer to OFs as 
regular generalized functions. The new objects (functions) that cannot be identified as OFs are called 
singular generalized functions (Think of the GF in Example 3 now.). Figure 3 shows how one should 
think of the space of GFs. In applications, several other test function spaces are used but the idea behind 
the construction of the associated GF spaces is exactly the same as we have presented above. At this 
stage, it is important to know that there are many singular GFs that are absolutely essential in applica- 
tions 7 . As a matter of fact, the introduction of the singular GFs into mathematics and the development 
of a rigorous operational algebra to work with them have provided a very powerful analytical tool to 
physicists and engineers. 



Singular GFs 


(J) g D 


G[(|>] =jg(x) <j)(x) dx 
g(x) symbolic function 

Real or complex 
numbers 


Figure 3- The space D' of GFs consisting of regular GFs (OFs) and singular GF that are new 


The readers should be aware that the definition of GFs given above, like the construction of real num- 
bers by Cauchy convergent sequences in analysis, is of enormous conceptual value. We will have to use 
this definition to obtain the properties of GFs and to extend the useful operations on OFs to all GFs. But 
once we have achieved these goals, we want to get away from dealing with functionals as functions and 
work with GFs with the same notations and ease as OFs. Like the man in the street who never needs to 
know about convergent sequences to work with numbers, we want to use GFs without being burdened 
by the conceptual framework of the construction of GFs. We can actually do this but first let us look at 
an example of a GF. 

Example 4- Here is a continuous linear functional defining a GF: 
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G[(f>] = 



e 2 x sin 3 x 0 (jc) dx - 5 0 (0) 


(peD 


( 8 ) 


You should have little problem to prove this fact. Note that, in practice, we always give the rule (the 
definition of the functional) only and we do not attempt to construct the table of the functional values 
because it is technically impossible. But we can use the rule to make part of the table if we want to! 

We will now dispose of the functional notation by the introduction of symbolic functions. This is a 
brilliant idea. It will allow us to work with singular generalized functions just like the OFs. We start 
with the Dirac functional of Example 3 which is a singular GF. We introduce a symbolic function 6(x) 
which is supposed to hold the information (or the memory) about the Dirac functional when we are 
doing algebraic manipulations until it appears as an integrand of an integral in a product with a test 
function. We then use the following relation agreed upon by convention : 



0 (x) 8(x)dx = 8[(f>] = 0 (0) 


(peD 


( 9 ) 


It is very important to recognize that ‘the integral’ on the left of this equation is not a Riemann or a 
Lebesgue integral because 6(x) is not an OF. You should think of the whole integral symbol like a 
Chinese character standing for d[0] = 0 (0) which is now completely meaningful. The symbolic function 
5{x), known universally as the Dirac delta function, can be used in algebraic manipulations like an OF 
until it appears under an integral sign like that on the left of eq. (9). You should, however, learn some 
more about the properties of singular GFs before you can comfortably or boldly manipulate them algebra- 
ically. We will give you what you need to know later but first let us look again at Example 4 using 
symbolic function notation. 

Example 5- The GF in Example 4 can be written symbolically as g{x) = e 2x sin 3 x - 5 8(x) . This GF is 
the sum of a regular and a singular GF. 


2.2 - Multidimensional GFs 


The above definition of GFs can easily extended to multidimensions. The test function space D is 
now n-dimensional and consists of all infinitely differentiable functions with bounded support. An 
example of such a function can be constructed from the function in eq. (4). Let x = (x\, X 2 , ...., x n ), 

and | x\ = (a'j + x\ + ... + x„) . Then the we can show that for a > 0, the following function is in space 

D: 


0 (x; a) = 


e \x \ 2 l{\x\ 2 -a 2 ) 

0 


| x | < a 
x I 2s a 


( 10 ) 


From this function and any continuous function, we can generate another test function by convolution 
over a finite region of space. Therefore, the space of test functions D has an uncountably many func- 
tions in it. We then define the multidimensional (Schwartz) GF space D' as the space of continuous 
linear functionals on D. This definition is appropriate for discovering the operational properties of GFs 
that give us a powerful tool for applications. In practice, however, we will work with symbolic functions 
when we deal with singular generalized functions. Therefore, we will be always working with general- 
ized functions as if we are working with OFs. Intuitively, this is of great help to people interested in 
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applications. The situation is similar to the use of identical notations for real and complex numbers as 
much as possible when we are solving a problem in complex plane. 

From our point of view, the multidimensional generalized functions are much more important than GFs 
of one variable in studying partial differential equations and in particular wave propagation problems. It 
is important to recognize that to study the multidimensional GFs, one needs to have a working knowl- 

in 

edge of differential geometry of curves and surfaces . 

As in the case of one variable, the most important singular GFs in multidimensions are the Dirac delta 
function and its derivatives. Of interest to us is the delta function with the support on a surface which 
can be in motion. This function has no one dimensional analogue. We start with a simple example and 
then discus the Dirac delta function with the support on a surface. 

Example 6- The Dirac delta function <5 (jc) with support at the origin has the following property 



0( x) 8(x)dx - 0 (0) 


(peD 


( 11 ) 


where V is an arbitrary volume which includes the origin. 

Example 7- Referring to Figure 4, let S : fix) - 0 be a smooth surface defined in such a way that the 
unit outward normal n = V/. This can always be done because if n ± V /, then define this surface by 
the new implicit function fix) / 1 v / | = 0 which has the desired property. 


f>0 n 



= Vf 

nit normal 


Figure 4- The description of the surface S defined implicitly by fix) = 0 with unit outward normal 

n = v/. 

We want to interpret the important GF 5{f) often appearing in applications. We define a Gaussian 
curvilinear coordinates (w 1 , ir ) over the surface S and define the third coordinate w 3 = / along the local 
normal to S to all points in the space in the vicinity of S . Note that we have 


d X- 


yj 8(2) («) d u 


U = (V , U 2 , M 3 ) 


( 12 ) 


where gg) (u) is the determinant of the first fundamental form of the surface 
u 3 = fix) = constant 29- 31 . We thus have the following very useful result 
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J' > 6(f) 0 (x) dx = J' i S(u 3 )(f) (x(u)) yj g( 2 ) (u) du l d u 2 d u 3 = 

(x(u l , u 2 , Oj) J g( 2 ) (u l , u 2 , 0) du l d u 2 = I f(x)dS 

Jf = o 


(13) 


You can see the importance of differential geometry here. Note that we have used the definition of the 
Dirac delta function in one variable when we integrate with respect to variable if after the second 
equality sign. For a more geometric derivation of this result see the papers by Farassat 29, 3 1 . 

It is important to remember that one cannot define a GF at a single point but over a finite interval. This 
does not cause any difficulties in applications. We can talk about equality of two GFs over an interval. 
We say that two GFs F[0] and G\d)\ are equal over an interval 1 if for all d> e D with support fully 
within I , we have F[(p] = G\d)\ . Remember, a functional is a rule that gives us a number for any function 
in D . Therefore, we are talking about equality of two numbers obtained from two rules. From this defini- 
tion we can say that 6 (x) - 0 (the zero function and not the real number zero) over the interval 
(-oo, -s\ IJ [e, oo) for any e > 0. 

To learn more about the multidimensional GFs we recommend the books by Gelfand and Shilov 7 , 
Stakgold 8 , Vladimirov 9 , Kanwal 11 , Jones 12 , Fenyo and Frey 19 , and the two NASA papers by Faras- 
sat 29,31 . For more advanced treatment as well some very important recent advances, see the books by 
Hormander 5 and Taylor 6 . 


2.3- Operations on GFs 

Now that we have extended the definition of ODs to GFs, we must extend the definitions of operations 
on OFs to GFs. The general approach is as follows. Write the definition of the operation in the language 
of linear functionals and then use it for all GFs if it makes any sense (see the definition of generalized 
Fourier Transform below where we have to change the test function space). In some cases the extension 
of the definition of an operation gives us the opportunity to solve new problems or solve old problems in 
much simpler ways. We single out generalized differentiation here as an example of an operation of 
considerable importance in applied mathematics. 

In the following we assume that all test functions are in space D . 


2.3.1- Simple Operations 

1- Addition of two GFs- If f(x) and g(x) are two OFs described by functionals F\d)\ and G\d )\ , d> e D 
then using the linear property of integrals, we can write 



+ g)(pdx 


J f <pdx+ J~, 


f (pdx + gfdx= F[0] + G[0] 


(14) 


Basically, we know the left side of the first equality is the functional representation of the sum of the 
OFs f(x) and g(x) and is given by the right side of the last equality sign. The above result now makes 
sense for all GFs, i.e., we define 


(F + G)[f] = F[f] + G[0] 


(15) 
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Note that (F + G) is a single functional which stands for the sum of the two GFs. 

Written in symbolic function notation, there is no difference between the notation for OFs and GFs for 
the sum of two functions. For example, for the sum of an OF fix) and the Dirac delta function we can 
simply write fix) + Six). 

2- Multiplication of two GFs- This operation cannot be extended to all of Schwartz GFs. In fact it is 
the aim of the new nonlinear GF theory to be able to define multiplication for all GFs 22,23 . Any two 
OFs can be multiplied. We can define multiplication of any GF with any C°° function. Again let fix) be 
any OF and let aix) e C °° , then 


/- 


-f 


aix) fix) fix)dx = I fix) iaix) 0(a)) dx = F[a 0] 


(16) 


We know that the left side of the first equality sign is the functional representation of the OF aix) fix) . 
Since the function afeD, the expression on the right of the third equality sign is a functional on the test 
function space D. The rule of multiplication of a generalized function F\d)\ and a C°° function aix) is 
given by the rule 


aF[f] = F[a 0] 


(17) 


Note that in this equation, the symbol aF stands for a single functional, i.e., a rule described by the right 
side of eq. (17). Again in symbolic notation we simply write aix) fix) but now we have a rule to inter- 
pret such products when fix) is a singular GF. We give an example here. 

Example 8- What is aix) Six) where a e C°° ? Let 0 e D, then 




a f>[0] = I aix) 0(A) dx = 5 [a 0] = a(0) 0 (0) = I a(0) 0(A) dx 


-!■ 


(18) 


We can write this result in symbolic function notation as 

aix) Six) = a(0) 0 (jc) (19) 


This useful result has been known for many years to physicists. 


3- Translation of a GF- Let fix) be an OF, and let us define E/ t fix) = fix - h ) , i.e., translating the 
function to the right for h > 0 . In functional notation, we have 




-f 


Eh fix) fix) dx- fix - h) fix)dx = fiy)fiy + h)dy = F\E_ h f\ 




( 20 ) 


Again the left side of the first equality sign is the functional generated by E/, fix) for which we use the 
symbol E/,F[f]. Since for 0 e D , we have 0 e D . Therefore, we use the following result for the 
definition of the translation of any GF: 


E h F[f\ = F[E_ h 0] 


( 21 ) 


Once again we have obtained a nontrivial result. In symbolic function notation we use the same notation 
for GFs as for OFs but occasionally we have to resort to eq. (21) to interpret the exact meaning of the 
translated GF such as in the following example. 
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Example 9- What is 6(x - a) where a is a real constant? The functional representation of this GF is 
obviously E a 8[f ] . Using the above rule, we have 



6(x - a) (p(x)dx = E a 0[0] = 6[E_ a 0] = 5[(p(x + a)] = 0(a) 


( 22 ) 


Again this result has been known by physicists for a long time. 

We can extend the definition of other simple operations on OFs to GFs such as the expansion and con- 
traction of the scale of the independent variable. For example, one can rigorously show that for a F 0 we 
have the following useful result 


1 

8{a x) = - — - 0(A) 
a 


(23) 


2.3.2- More Advanced Operations 

We will extend the definition of two operations on OFs to GFs here. These are Fourier transformation 
and differentiation operations. 

1- Fourier Transform of GFs- We work in one dimension here. Let f{x) be a function defined on IR 
which has Fourier transform define by the relation 



fix) e 2nixy dx 


(24) 


If the functional F[0] is identified with fix ) , then its Fourier transform must be identified with f(y) as 
follows 


/~*oo /~*oo 

F[(p]= I f(y)<f>(y)dy= I f(y)$(y)dy = F[(f>] (25) 

where 0 is the Fourier transform of 0. We have used Parseval’s theorem after the second equality sign. 
Now the functional F[0] is only defined if 0 e D whenever 0 e D . We can show that, in general, this is 
not so and 0 can lie outside the space D. We need to define a test function space S on IR consisting of 
C°° functions which go to zero at infinity faster than \x~ n \ for any n . This can be done and in fact we 
can show that D c S. The space of tempered GFs S' i s the space of all continuous linear functionals on 
S. An important relation to remember is D c S c S' c D' . All the GFs in the space S' have generalized 
Fourier transform given by the relation F[0] = F[0] where now feS. We will say no more on this 
interesting subject here. See the references at the end of this paper. We strongly recommend the book by 
Strichartz 10 for this subject. 

2- Differentiation of GFs- This is the most important operation on GFs because of its many applica- 
tions in mathematics. Many of the problems of wave propagation need the use of this concept. Let us 
start with a C 1 function fix) which we identify with the functional F[0] . The functional representation 
of fix) can be manipulated as follows 
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F'\4>\ 



'(.*) 0 (x) dx 



(jc) dx = 


-Fm 


(26) 


Note that we have integrated by parts after the second equality sign and we have used the fact that the 
support of 0 is compact to drop the terms involving the limits of the integral. Since we know that the 
function 0' is in D, F[f'] is a proper functional on D and we can use the following relation for the 
definition of the generalized derivative (GD) of a GF 


F'm = -Fm 


(27) 


From this we can define the nth generalized derivative of a GF 

F {n) [(t>] = (-1 ) n F[(p {n) ] 


(28) 


This significant result states that GFs have GDs of all orders. When we work with symbolic functions 
we often utilize a bar over a derivative operator to designate generalized differentiation if the 
operation can be confused with ordinary differentiation. See eqs. (30) and (34) below. It is understood 
that the differentiation of singular GFs can only be GD and in that case we will not use this notation. 

Example 10- GD of the Heaviside function h{x ) = 0 lor x < 0 and h(x) = 1 for jc >0 is found as follows 

FI' [(f) ] = -H[(f)'] = - (f)' (x) dx = (f) (0) (29) 

Jo 


Note that we have used the fact that the test function has a bounded support and, thus, we have 
0(oo) = 0 . Therefore, we have 


h'(x) = 0(x) 


(30) 


That is, the GD of the Heaviside function is the Dirac delta function . Note that h'(x) = 0 but h\x) + 0. 
We will return to this point later. This result shows that the GD of an OF can be a singular GF! 

Example 11- Generalized derivative of the Dirac delta function is 


6'[(f>\ 



(x) (p (x) dx 


-sm 


-0'(O) 


(31) 


Similarly, we have 


0 W [0] 



6 {n \x)(p{x) dx= (- 1)” 0 ( ")(O) 


(32) 


Note that in eqs. (31) and (32), the integrals are meaningless and stand for the functional rules to the left 
of the first equality sign in the respective equations. 

Example 12- Let the OF f(x) be differentiable with a single jump of A / = f(x o+) - /'Go- ) at the point 
.*0 as shown in Figure 5. 
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Figure 5- A differentiable function with a single jump discontinuity at xo 

We find the GD of this function as follows. Let F\(f>\ be the functional representation of this function. 
Then its GD can be manipulated as follows 


F'[4>\ = -F[(f) 




f{x) (f)'(x)dx 


nx 0 - 

— OO 


fix) fix) dx 


X) 


fix) fix) dx - 


r»*o- 

J ? 

V/ — OO 

r*oo 

J r 

— OO 


(x) fix) dx + 


r*oo 

J r 

«T*o+ 


(*) fx) dx+ (fix 0+) - fix o-)) f(x 0 ) = 


(33) 


(x) fix) dx + A / fix 0 ) 


r*00 

= (/' 

— OO 


ix) + A f - xq)) f(x)dx 


Interpreted in symbolic function notation, we have obtained the following important result 

f\x) = f(x)+ Afdix-xo) 


(34) 


Again we have shown that the GD of an OF can be a singular GF. Furthermore, this result shows the 
need for a notation for generalized differentiation which we have utilized here, i.e., the overbar on the 
left side. 

Example 13- Referring to Figure 4, let the function k(x) be a differentiable function with a jump discon- 
tinuity across the surface f = 0. We define the jump of the function k(x) across the surface as follows: 

Ak(u l , m 2 ) = k(u l , u 2 , 0+) -k(u l , u 2 , 0_) = kix \f=o+) - k ix l/=o-) (35) 

where u = (w 1 , ir , id ) is the curvilinear coordinate system we used in Example 7. Note that the func- 
tion k ix) is discontinuous in variable ir = f only. This is precisely why we selected such a coordinate 
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system in the vicinity of the surface of discontinuity. Now we are going to use eq. (34) to derive the 
following important result: 


dk(x («)) 
dxi 

dui dk 


dx; duJ 


, du 3 <•( dk , df c- r dk , 
r + A k — — 8 (m 3 ) = — \- A k — — 8 (/) = — A k n, 8 (/) 


(36) 


dxi 


dx; 


dx; 


dxi 


We have used the summation convention on repeated index in the above equation. We can write the 
above result in vector notation as follows: 


V k (x) = V k (jc) + Akn 8 (f) 


(Note n = v/here.) 


(37) 


We give here two more important relations which can be remembered easily. Let k (x) be a vector field 
with a jump discontinuity across the surface / = 0 with the jump defined similar to eq. (35). Then the 
following results hold: 

(Note n = V / here.) (38) 


(Note n = v / here.) (39) 


V • k (x) = V • k (x) + A k ■ nd(f) 


vxk(x) = V x k (x) + n x A k 8 (/) 


With what we have obtained so far we can solve a lot of interesting problems. Before we give some 
examples, let us say a little more about why one should use generalized derivative. 


2.4- Why We Use Generalized Derivatives 

The most important fact to remember about GDs is that in GF theory generalized differentiation is a 
continuous operation. In particular, generalized differentiation and integration are inverse operations. 
This means that the integration of GD of a GF recovers the function faithfully with all of its jumps. 
Examine the following 

Example 14- Let h(x) be the Heaviside function defined in Example 10. Then we can easily show that 

r*x r*x 

I h'(y)dy = h(x ) but I h ' (y) dy = constant + h (x) (40) 

J —co J —oo 

Therefore, any problem involving unknowns with discontinuities should be set up in GF space and 
all derivatives must be treated as GDs from the start. Such discontinuities are either artificially 
introduced, e.g., as in FW-H equation, or are genuine, e.g., a shock wave in the flow field. The state- 
ment following eq. (40) has profound implications for applications some of which will mentioned below: 

1- All local conservation laws obtained from control volume analysis are valid if all the derivatives 
are considered GDs. Examples of such conservation laws are the mass continuity and the Navier- 
Stokes equations. 

2- The jump conditions across a natural discontinuity in the flow are inherent in the conservation 
laws. These jump conditions are obtained easily by applying the definition of GD to the conservation 
laws. See Example 15 below. 
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3- Any problem involving derivatives in which singular GFs appear must be set up in GF space. 

One must then carefully distinguish between the GDs and the ordinary derivatives of any function 
appearing in the algebra. An example of such a problem is finding the Green’s function of an ODE or a 
PDE. 


4- Generalized differentiation operator commutes with all common limit operations. This is the 
most pleasant and useful property of generalized differentiation. The following operations are allowed 
provided that the limits of the integration is independent of the variable of differentiation: 


d 


dx 


~ I q(x,y)dy= I 
i J £1 J n 


dq(x,y ) d dq n (x,y) 

n dy, —— lim q n (x, y) = lim 

OXi OXi n-* OO >i->00 OXi 


(41) 


See the NASA paper by Farassat 3 1 for more examples. 

Example 15- Let us find one of the jump conditions across an unsteady shock wave which is described 
by the surface fix, t) - 0. We define the jumps in fluid parameters as in eq. (35). The mass continuity 
equation is 


dp_ 

dt 


+ V • (p m) = 0 


(42) 


Here the fluid density is p and the fluid velocity is u. First we know this conservation law was derived 
from a control volume analysis. Therefore, it is valid if we treat all ordinary derivatives as GDs, i.e., we 
have 


dp _ 

— + V • (p k) = 0 
Ot 

We now use the definition of GD in this equation as follows 

dp _ 

— - + V • (p U) = 

Ot 

■— + A p ~~ 6 (/) + V • (p u) + A (p u) ■ n 5 (/) = A[p (u n - v n )] 8 (/) = 0 
ot ot 


(43) 


(44) 


where we have used the notations u n = u n and v n = -df/dt is the local normal velocity of the shock. 
We have used eq. (42) in the above result. Therefore, one jump condition across an unsteady shock 
wave is 


A[p (u n - vj] = 0 or A (p u n ) - v n A p = 0 (45) 

Note that we derived this condition without using the pill box analysis employed in the classical deriva- 
tion of the result. 

5- We can use the Green’s function technique to find the discontinuous solutions of an ODE or a 
PDE provided that we formulate the problem in GF space and use generalized differentiation 
everywhere. This is a very powerful result and its impact on applications is enormous. We will not give 
any example here but you can see some good examples in the next section. 

We conclude this section by stating the fundamental theorem that characterizes GFs in D ' : 
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6- GFs in D' are GDs of finite orders of continuous functions 7 . This is a very profound result 
obtained by Laurent Schwartz. As a simple example we note that the Dirac delta function is the general- 
ized second derivative of the continuous function that is zero on the negative axis and is equal to x on 
the positive axis. 

2.5- A Note on Test Function Spaces 

We have developed above the GF theory using the test function space D. Although the GF spaces D' 
and S' are very useful in applications, we often have to use test functions other than the spaces D and S 
which are not C°° functions. But the construction of GFs in the space D' can be used as a model to 
get other kinds of GFs. For example if we take a space of test functions consisting of C n functions, 
then we can define GFs on this space as continuous linear functionals which only have GDs up to the 
order n. We must, of course, modify our definition of the continuity of a linear functional appropriately. 
To define GFs suitable to treat ODEs and PDEs, say a differential operator L (a differential equation 
plus linear homogeneous boundary conditions (BCs)), the test functions must have continuous deriva- 
tives of some order (related to the order of highest derivative of L) and should satisfy the adjoint BCs. 
Everything we have learned about GFs here, particularly the operations on GFs, apply as they are. For 
this reason we do not specifically identify the test function space unless such identification will shed 
some light on the problem at hand. 

3- Some Applications of Generalized Differentiation 

We will give two important applications of generalized differentiation which are very useful in solving 
wave propagation problems. These are: 

1- The Green’s function of a second order ODE, and 

2- The imbedding of a problem in another problem whose Green’s function is known. 

You should note that in these applications we have to utilize almost everything we presented on GFs in 
the previous section. 

3.1- The Green’s Function of a Second Order ODE 

In working with wave propagation problems in time or frequency domains, we often have to obtain the 
Green’s function of a second order ODE in intermediate steps. Since we are now able to explain some of 
the subtleties of the derivation of this Green’s function from the point of GF theory, we will present a 
discussion of it here. We hope that you will agree with us that GF theory is the right tool to use here. 

Consider the following second order ODE with two linear homogeneous boundary conditions (BCs): 

l u — a(x) u"(x ) + b(x ) u'(x ) + c(x ) u{x) = k{x) x e [0, 1] (46) 

BCj [u] = 0, BC 2 [u] = 0 (47) 

The Green’s function g(x, y) for this problem is formulated classically as follows 8 (see volume 1 of this 
reference): 
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lx g(x, y) = a(x ) — + b(x) ~ + c(x) g(x, y) = 6 (x- y), x, y e [0, 1] (48) 

OX A ox 

BCi[g(x, y)] = 0, BC 2 [g(x, y)] = 0 (in variable jc only) (49) 

Now we see that we have the Dirac delta function on the right of eq. (48) which is a singular GF. So the 
partial derivatives in this equation are most likely GDs (we know the answer but we are presenting what 
most books for engineers and physicists give). But how did we suddenly get into the GF space? Well, 
we did not really set up our problem correctly. First from the theory of ODEs we must know something 
about the nature of the solution of the BC problem described by eqs.(46) and (47). Let us assume that 
we know that ueC 2 . This can be proved a posteriori when the solution is at hand. In that case, in the 
domain of the definition of this function, since a C 2 function cannot have any jumps or sharp corners, 
we have 

~lu = a(x ) u"(x) + b(x ) ~u' (x) + c(x) u(x) = a(x) u"(x) + b(x) u'(x ) + c(x) u(x) = l u (50) 
The Green’s function is used to find the unknown function from the following relation: 

u(x)= I k(y) g(x, y)dy (51) 

Jo 


Let us discover further properties of the Green’s function from this relation by applying the linear differ- 
ential operator (the ODE plus the BCs) to both sides of the equality. First we get 


/ u(x) = / u{x) 


-T 


k(y) g(x, y)dy 


-f 


k(y) l x g(x, y)dy = k{x) 


(52) 


Note that we freely took the generalized derivatives in l x inside the integral sign without concern 
because this step is allowed in GF theory. From this we immediately can see that 


lx g(x, y) = a(x) y -- + b(x) d§ ^ C, - y - > - + c(x) g(x, y) = 6 (x - y) 

OX L ox 


(53) 


Therefore, we established rigorously that all the derivatives in eq. (48) are GDs. We will get back to 
the implications of this fact soon. Next we apply the two BCs to both sides of eq. (51) remembering that 
these are linear and homogeneous: 

BC X [u(x)] = BC a f k(y) g(x, y)dy = 

Jo 

where the symbol BC X stands for either of the two BCs but applying to variable x only. From the above 
result, because the function k(y) is arbitrary, we conclude that the Green’s function in variable x 
satisfies both BCs: 


/ 


k(y)BC x [g(x, y)]dy = 0 


(54) 


BCi ?x [g(jc, y)] and BC 2 ,.Jg(x, j)] 


(55) 
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Now let us see what eq. (53) implies about the Green’s function. Obviously this function changes its 
character at x = y . So let us assume that it is described as follows by two C 2 functions with possible 
jumps in their values and their derivatives at x = y: 


g( ^ ) = W,y) 


x < y 
x > y 


(56) 


First of all since the Dirac delta function S (x- y) is identically zero in the two open intervals on both 
sides of x = y, we must have: 


hg\(x, y) = l x g 2 (x, y) = 0 


(57) 


This means that the functions g\(x, y ) and g 2 (x, y ) in variable x are solutions of the homogeneous ODE 
l x g(x, y) = 0. Let us take the GD of the first and second order of the Green’s function described by eq. 
(56) 


dg(x,y) dg(x,y ) 


dx 


dx 


+ [g 2 (j+, y) -gi(y~, y)]S{x-y) 


d g(x,y ) d 2 g(x, y) 


dx 2 


dx 2 


+ 


dg2(y+,y) dgi(y-,y) 


dx 


dx 


(58) 


(59) 


S (* - y) + [g2(j+ , y) - gi (y- , j)] S' (x - y) 


We must substitute these two relations in eq. (53) and equate both sides term by term. Since there is no 
derivative of the Dirac delta function on the right of eq. (53), we immediately conclude that 


g2(y+, y)-gi(y~, y) = 0 


(60) 


This means that the Green’s function is continuous at x = y. Furthermore, we get 

(assuming that a(y) + 0 fory e [0, 1]) (61) 


r dg 2 (y+,y) dgi(y-,y)~ 

1 

dx dx 

a(y) 


This means that the derivative of the Green’s function with respect to x has a jump described by the 
above equation at x = y. We can show that eqs. (55), (57), (60) and (61) uniquely define the Green’s 
function. Note that only by using GF theory we have been able to get the above results rigorously and 
satisfactorily. 

Example 16- Let us take / u{x) = u" for xe [0, 1] with the following BCs: 

m(0)-2m'( 0) = 0 and m(1) + m'(1) = 0 

The above procedure will give the following Green’s function for this problem: 


(62) 


g(x, y) = 


( (v-2)(.r + 2) 
4 

(v + 2) (x - 2) 
4 


x < y 
x > y 


(63) 
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This Green’s function satisfies the conditions derived above. Note that the differential operator in this 
example is selfadjoint and, therefore, the symmetry of the Green’s function in variables x and y is 
expected 8 . 

3.2- The Imbedding of a Problem in Another Problem 

We got interested in this problem through our work on the Ffowcs Williams-Hawkings (FW-H) equa- 
tion 2 . The situation is described in Figure 6 below. 



this data surface 

Figure 6- The problem of noise radiation from a moving body (a helicopter rotor blade here) as 
formulated by Ffowcs Williams and Hawkings 

We are given an inhomogeneous wave equation in the exterior of a moving surface described implicitly 
by the function /( x, t) - 0 . This surface will be called the data surface here. Ffowcs Williams and 
Hawkings searched for a radiation formula in the exterior of the data surface. Essentially, these authors 
used an imbedding method by making their wave equation valid in the entire unbounded space^ . This 
imbedding allowed the use of the free space Green’s function for getting the radiation formula. They 
assumed that the fluid inside the data surface is always at rest with the condition of the undisturbed 
medium. Since an artificial discontinuity has been introduced in the medium, one must work with GFs. 
This imbedding technique provides us with an elegant method of deriving a radiation formula. This 
technique can be applied to a whole range of other problems involving ODEs and PDEs particularly 
those related to wave propagation. We will discuss some of these problems here. First, however, we 
present the imbedding technique. 
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3.2.1- The Imbedding Technique 

We will consider the problem of radiation from a moving and deformable surface f(x, t) = 0 , assuming 
that v / = n, where n is the unit outward normal. We have shown the moving surface at two times in 
space in Figure 7. 



n= Vf unit normal 


Figure 7- A moving deformable surface described by f{x, r) = 0 at two different times t\ and t 2 in 
space 

We want to solve the following exterior acoustic radiation problem: 

1 d 2 (f 7 

V z (fi = 

c 2 dt 2 (64) 

□ (p = Q(x, t) ( x in exterior of / = 0, c = constant) 


We are assuming that we have all the data about the function ip on the moving surface. These data may 
be overdetermined but right now we are not concerned about this fact. We want to extend the domain of 
the problem to the interior of the surface so that we can use the simple free space Green’s function of 
the wave equation. In other words, we want to imbed our problem in another problem with larger 
domain in space (this can be done also for the dimension of time). This imbedding will require the 
powerful machinery of GFs. Here is how we do the imbedding. We start by extending the definition of 
the two functions in eq. (64) as follows: 


<P (x, t) = 


V (x, t) 
0 


/> 0 

/< O’ 


and 


Q (x, t) = 


f Q (x, t) 

10 


/> 0 
/< 0 


(65) 


The first thing we notice is that with this extension, the new functions satisfy the following wave equa- 
tion: 


□ 2 <p = Q (x, t) 


( 66 ) 
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Note that the unknown function p is discontinuous. To find this function, we must set the problem in 
GF space. This simply means that we find what the wave operator looks like when all the ordinary 
derivatives are written as GDs as follows: 


dip 

dt 


dp 

dt 


+ <p 


df_ 

dt 


6(f) = 


dp 

dt 


<P V n 6(f) 


(67) 


where v n = - df/dt is local normal velocity of the surface which is unambiguously defined even for a 
deformable surface. Also we use the notation <p = <p(f = 0+) when this function multiplies the Dirac 
delta function 5 (/) as in the last term on the right of eq. (67). Next we find the following relations: 


d p 
dt 2 


d 2 p 
dt 2 


dip 

dt 


V„ 6(f) 


d 

dt 


[ip v n 6(f)\ 


( 68 ) 


V 2 p = v 2 p + n-Vip 6(f) + v-[ipn 6(f)] 


(69) 


where we have again used the notation that any function of p multiplying the Dirac delta function is 
evaluated at / = 0+ . Let us define the local normal Mach number of the surface as M n = v„/c, and 
define additionally the following two symbols tp t - dp / dt and <p n - n ■ v <p. Then the above two rela- 
tion will give us: 

□ 2 p = D 2 p - — ip t M n 6(f) - — — [ip M n 6(f)] -ip 6(f) - V • [ip n 6(f)] (70) 

c c dt 


Or, after using eq. (66), we obtain 


□ z P = ~ M n 6(f) - — —[ipM n 6(f)] - ip 6(f) -v-[ipn 6(f)] 
c c dt 


1=:2 


(71) 


Note that we do need to use a bar over the wave operator on the left side since from the right side it is 
clear that all derivatives involved are GDs. Now we can use the Green’s function of the wave equation 
in the unbounded space, the so-called free-space Green’s function, to find the unknown function 
p{x , t) everywhere in space. The result is the Kirchhoff formula for moving surfaces 31 ' 32 . The 
classical derivation of this formula by Morgans 33 is very complicated and in fact until Farassat and 
Myers rederived it by the modern method presented here, there were some doubts expressed about the 
correctness of Morgans’ result 2,34 . 

We will not derive the Kirchhoff formula for moving surfaces here since the two references by Farassat 
and Myers 32 and Farassat 31 are quite comprehensive alleviating the need for further elaboration. How- 
ever, we next derive some classical results for stationary surfaces by the imbedding method which are 
surprisingly simpler than the classical methods of their derivation. 

Before we end this discussion we want to mention an important question that comes to mind. Answering 
this question may lead to further applications of the Green’s function of the unbounded space. To imbed 
our exterior radiation problem in IR 3 we have imposed that inside the surface both p and Q be zero. 
Can the functions p and Q be assumed to take other values and are there advantages in doing so? 
We will find the answer to this question using an ODE to simplify the analysis. 
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3.2.2- Derivation of Some Classical Results Using the Imbedding Technique 

1- The Green’s Theorem for Laplace Equation in the Exterior Domain of a Surface- The problem 
is described in Figure 8. 



Obs. 


Figure 8- Problem description for deriving the Green’s theorem for Laplace equation, / > 0 out- 
side the surface and n - v / on the surface 

We imbed this problem in IR 3 by extending the definition of tp to the interior of the surface as in eq. 
(65). We find generalized Laplace equation next: 

V 2 <£ = V 2 p + n-Vipd(f) + V-[ipnd(f)] = tp n 6(f) + V ■ [ip n 6(f)] (xeR 3 ) (72) 


The definition of all symbols are identical to those in Section 3.2.1. We have used the fact that v 2 <p = 0 
in IR 3 . The Green’s function of this problem in the unbounded space is -1/(4 nr), where r = | x — y \ . 
Here x and y are the observer and the source (field) variables, respectively. The solution of the above 
equation will give us the answer we are seeking as follows: 


<p(y)n 


4n(p(x) = - f y - 8 (f)dy -v • f ------ - - 5 (f)dy = 

Jr 3 r *7ir 3 


-/ 


r 3 r 

<Pn (30 

r 3 r 


r 3 r 

1 


$ (f)dy- if(y)n-v x 1 8{f)dy = 


(73) 


U- 


^ n (y) tp cos 6 


+ 


dS 


We have used the result of Example 7, i.e., eq. (13), twice in the last step. The angle 6 is between the 
local outward normal to the surface and r = x- y, see Figure 8. This is the well-known Green’s theo- 
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rem for the Laplace equation. Note that a bonus of the imbedding method is that we know that by con- 
struction this equation will give <p = 0 inside the surface, a fact not at all obvious in classical derivation 
and must be established separately. 

2- The Classical Kirchhoff Formula (The Exterior Problem)- We will now give the derivation of this 
well-known result based on eq. (71), i.e. using the imbedding method. Referring to Figure 9, we imbed 
the problem in R 3 and use the fact that M n - 0 in eq. (71) to get: 

n 2 (p = -<p n 6(f)~ v-[<pn6(f)\ (74) 


X 



(Kirchhoff) data surface 

Figure 9- Problem description for deriving the classical Kirchhoff formula, / > 0 outside the 
surface and n = V / on the surface 

The Green’s function of the wave equation in the unbounded three dimensional space is: 


G(x, f; y, r) = 


0 T > t 

6(r-t + r/c)/4nr r^t 


(75) 


where r = | x - y | . Here (x, t) and (y, r) are the observer and the source space-time variables, respec- 
tively. We usually use the following symbol: g = t — t + r / c. 

We are first interested in the solution of the equation: 

□ 2 if/ = q (x, t) 6(f) (76) 


The solution of eq. (76) is 8,9,29 : 

4 ji if/(x, t) = J' q (j, r)S(f) dydT = J~ 


q(y,t-r/c) 


dS 


f = o 


(77) 


The limits of the integral in this equation are given below: 
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/ /~*l r* r*t r* CX) /~*oo s*oo 

dydr- I I dydr- I I I I dy\ dy 2 dyi dT 

%J — OO^tR ' 1 —oo —oo —oo —oo 


(78) 


-oo *7|R 

Therefore, the solution of eq. (74) is: 

(y, t-r/c) 

>f = o r 


4 7T (f (x 


, 0 = - f 

Jf=( 


dS-V , 


f 

J/ = o 


ntpiy, t-r/c ) 


7S = 


<p(y, t-r/c) 


- f <Pn ^ y,t r / - c) dS- f n-V x ( 
J/ = 0 r J f = 0 ' 


'/ = o 


(79) 


Taking the partial derivatives with respect to the variable x of the integrand of the second integral, we 
obtain the following result: 


4 71 (f (x, t) = 

c~ l cos 6 (ft (y, t-r/c)-(f n (y, t-r/c ) 

7 = o r 


l 


dS + 


f 

J/=o 


cos 6(f(y, t-r/c) (80) 
^ d 


The angle 6 is between the local outward normal to the surface and r = x - y, see Figure 9. This is the 
well-known classical Kirchhoff formula 34- 36 . Note that we know by construction that (p = 0 inside 
the surface. This fact is not obvious in the classical derivation and, as in the case of the classical deriva- 
tion of the Green’s theorem for the Laplace equation, must be established separately. Notice also that in 
our derivation we did not have to resort to the four-dimensional Green’s theorem for the wave equation. 
This is a great advantage of the imbedding method over the classical derivation. 

3- Radiation into Half Space- The Rayleigh Integral- We will derive two well-known results for 
radiation into half space by the imbedding method. Let the half space be defined by > 0 . We are, 
therefore, interested in solving the following problem: 

n 1 t p = Q(x,t) (x 3 >0) (81) 


We will say nothing about the BC now because two possible BCs emerge from the imbedding method. 
Let us imbed this problem in R 3 in two ways as follows: 


<P (x, t) = 

( <P (x, t) 
l ~(f (x, t) 

JC 3 > 0 
JC 3 < 0 ’ 


Q(x, 

(Q (x, t) 

; “I -Q(X, t) 

JC 3 > 0 
*3 < 0 

(Imbedding 1) 

and 




<P (x, t) = 

( <P (x, t ) 
l <f (x, t) 

*3 > 0 
*3 < 0 ’ 


Q(x, 

)_ iQ (x, t ) 

\Q(X, t) 

*3 > 0 
*3 < 0 

(Imbedding 2) 


(82) 


( 83 ) 
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Imbedding 1- We can show easily that on the plane .* 3=0 the jump 

A ip — 2 Lp (x [ X 2 , 0+ , t) = 2 q (x [ . a' 2 , t ) , and also A (dtp/ dx?,) = 0. Therefore, the extended function 
satisfies the wave equation: 


(xeU 3 ) 


□ 9 (p = Q (x, t) - 2 q Op x 2 , t) S'(x 3 ) 

Using the free space Green’s function of the wave equation, we get 

Q(y, t-r/c ) r f 2S'(y 3 )q(y h y 2 , t- r / c) 

/R 3 r 


(84) 


4 7T ^ (JC. 


,t)= f 

J IR 3 


L 


Q(y, t-r/c ) 


dy + 


f 

J IR 2 


dy- f 

J IR J 

d 


dy = 


- 


2q , (ji. J2, t-r/c) 


(85) 




"/3 = 0 


When we take the derivative in the integrand of the second integral, we obtain the following result: 


4 7T Cp (JC. 


,t)= f 

J IR 3 


Q(y, t-r/c) 


dy + 


L 


( c 1 cos0q t (y\.y 2 , t-r/c) cos 6 q (y L y 2 , t - r / c) 


( 86 ) 


+ 


dy\dy 2 


Here we have used q, -dq/dt and the angle 6 is between the radiation direction x-y and the 
x 3 - axis. Note that the function Q is defined in eq. (82). This result is the solution of the following 
problem in the half space x 3 > 0 : 

(*3 > 0) (87) 


□ 2 (fi = Q (X, t), if (xi t x 2 , 0 ,t) = q (x h x 2 , t) 


Imbedding 2- We can show that on the plane .*3=0, the jump A ip — 0 , and also 
A(d<p/dx 3 ) = 2d<p(xi'X2,0+, t)/dx 3 = 2q(xi'X 2 , t). Therefore, the extended function satisfies the 
wave equation: 

U 2 = Q(x, t)- q(x h x 2 , t)S(x 3 ) (jcefR 3 ) (88) 


Using the free space Green’s function of the wave equation, we get 

, ^ f Q(y, t-r/c) f 2d(y 3 )q(y L y 2 , t-r/c) 

4 n<p(x,t)= dy- 

Jr 3 r Jr 


dy = 


IR' 


f g(j> t-r/c) dy _ C r 2q(yi,y 2 , t-r/c) 
Jr 3 r J R 2 [ r 


(89) 


Jv 3 =0 


dy\dy 2 


Since the source in the last integral lies always on the V 3 = 0 plane, we can drop this subscript in inte- 
grand and we obtain: 
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4 7T (f (x, t) = 



Q(y, t-r/c) 
r 



2q(y u y 2 , t-r/c) 
r 


dy i dy 2 


(90) 


Note that the function Q is defined in eq. (83).This is the solution of the following problem in the half 
space JC 3 > 0 : 


(*3>0) (91) 


Equation (90) without the volume term is known as the Rayleigh integral 36 . 

4- The Imbedding Problem for a Second Order ODE- We consider the second order ODE 
1 u(x) = u"(x ) - k(x) for xe[a, b ] c [0, 1] with two linear and homogeneous BCs which we denote as 
BC'i [»] = 0 and BC^M = 0. Assuming that the solution exists, we want to find the solution by imbed- 
ding method using the Green’s function of the Example 16. The situation is shown in Figure 10. 


□ <P = Q (x, t ), 


dip 


(x h x 2 , 0, t) = q(x h x 2 , t ) 



0 a b 1 


Figure 10- Imbedding an ODE problem for the interval [a, b ] in the larger interval [0, 1] for the 
same differential equation 

Note that the ODE of Example 16 is the same as that the one on the smaller interval This is required in 
an imbedding method. Let us extend the unknown function u(x) to the larger interval as follows 


[ u(x) xe[a, b] ~ f k(x) xe[a,b] 

< , and k(x ) = < 

l0 xe [0, 1] \[a, b] \0 xs [0, 1] \[a, b] 


(92) 


See Figure 10 above. Notice that u satisfies the BCs of Example 16 given by eq. (62). This is required 
because if this condition is not satisfied, then we cannot use the Green’s function of Example 16 given 
by eq. (63). Note that we have / u = u" = f . We next find the following result valid on [0, 1] : 


1 u (x) = 

u\x) = u"(x) + u'(a) S (x — a) - u'{b ) 8{x — b) + u (a) 8\x -a) - u(b ) 6\x - b) = (93) 

k(x) + u'(a) 6 (x - a) - u'(b ) 8 (x- b) + u (a) 8'(x - a ) - u(b) 8'(x - b) 


Now let g(x, y ) be the Green’s function of Example 16 given by eq. (63), then we can find a representa- 
tion of the solution as follows: 
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u(x ) = 


f k(y) g(x, y) dy + u'(a) g(x, a) - u'(b) g(x, b) + u (a) — ^ — -- uib ) ^ ^ 

Jo 


dy 


dy (94) 


r h , , , . dg(x,a ) dg(x, b) 

I k(y) g(w y)dy + u (a) g(x, a) - u ( b ) g(x, b) + u (a) u{b ) 

Jci 


dy 


dy 


We are not finished yet. We can show that using eq. (94), w(a_) = 0 and u(a + ) = u(a ) give identical 
relations among the boundary values (BVs) of the imbedded problem. Similarly, we can show that 
u{b+) = 0 and u(b - ) = u{b) also give us another identical relations among the BVs of the imbedded 
problem. We also have two relations supplied by the BCs of the imbedded problem. These four relations 
will give us the four unknowns u(a), u(b), u'(a), and u'(b) which we can then substitute in eq. (94) to 
completely solve this problem. 

Looking at this procedure we note that all we have required of the behavior of the extended function 
u(x ) is that it should satisfy the BCs of the problem whose Green’s function is available in the larger 
domain. This means that any solution v(x) of Example 16 can be assumed in the extended region. 
Assuming that we have l v(x) = j(x ) , then eq. (92) must be modified as follows: 


u{x) = 


I u(x) 
lv(jc) 


xe [a, b] 
i£[0, 1] \[a, b] 


, and 



xe [a, b] 
xe [0, 1] \[a, b] 


(95) 


and eq. (93) becomes 

l u(x) = 

u\x) = u'(x) + [ u (a ) - v'(a)] 6 (x - a) - [ u'(b ) - v'(b )] 6 (x - b) + [u (a) - v(x)] 

S' (x - a) - [u ( b ) - v(b )] S'(x - b) = (96) 

£(x) + [ u'(a ) - v'(fl)] S (x - a) - [ u'(b ) - v'(b)\ S 

(x - b) + [u (a) - v(x)] S' (x - a ) - [u (b) - v{b)\ S'(x - b ) 


Again using the Green’s function of Example 16, we find the representation of the unknown function as: 


u(x)= I k(y)g(x,y)dy + [u'(a)-v'(a)]g(x,a)-[u'(b)-v'(b)] 
Jo 


g(x, b) + [u (a) - v(x)] 


dg(x, a) 
' ~dy 


[u (b) - v(b)] 


dg(x, b) 

dy 


(97) 


Two relations are obtained from u{a+) = u{a) and ii(b-) = u(b) among the BVs of the imbedded prob- 
lem. Together with the two BCs of the imbedded problem, we can then find u(a), u(b), u'(a), and u'{b) 
and substitute in eq. (97) to find the solution of the imbedded problem. 

The implication of this result for wave propagation problems can be substantial. It means that we do not 
have to assume that fluid inside the moving body is at rest with the conditions of the undisturbed 
medium in the derivation of the FW-H equation. We can take any solution of the wave equation 
□ 2 ip = q inside the surface and then apply the imbedding technique. This will allow us to use the free- 
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space Green’s function and perhaps not have some of the shortcomings of the FW-H equation such as 
the cancellation between the contributions of the surface and volume integrals. 


4- Concluding Remarks 

In this paper we have discussed some applications of GF theory principally to wave propagation prob- 
lems. We have first given a fairly complete tutorial on GF theory and then have utilized most of the 
properties of GFs in the examples and applications presented. We have emphasized generalized differen- 
tiation here which is a very important operation on GFs with many applications. These are some of our 
main conclusions: 

1- GF theory is a very powerful analytic tool in physics and acoustics with a wide range of applications, 

2- The derivation of many classical results in physics and acoustics can be enormously simplified by 
using GF theory, 

3- The power of the Green’s function technique is substantially increased if GF theory is utilized as in 
the imbedding technique, and 

4- It is possible to use functions other than the null function inside a moving surface in the imbedding 
problem used for obtaining the FW-H equation. Some of these choices may alleviate the shortcomings 
of the FW-H equation and still allow us to use the free-space Green’s function. 

Generalized function theory is a very powerful tool in mathematics that every physicist and mathemati- 
cally minded engineer should strive to learn. We have tried to impart to our readers the pleasures of 
learning and using GFs in their work. 
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